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Is the resonance D(2637) really a radial excitation? 
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Abstract 

We consider various possible identifications of the quantum numbers of the resonance D(2637) recently 
observed by DELPHI in the D*tttv channel. We argue that in spite of a good agreement of the measured 
mass with the quark-model prediction for the radial excitation, a total width as small as < 15 MeV is hardly 
compatible with its identification as a radial charm excitation. The J p — 2~,3~ orbitally excited mesons 
with such a mass could have widths of the observed order of magnitude. However in this case one would 
ON ' expect two neighbouring states with the mass difference of about 30-50 MeV corresponding to the nearly 

degenerate components of the heavy-meson multiplet with light-quanta angular momentum j=5/2 , and 
fl j' moreover, according to the quark-model predictions the mass of the orbital excitation should be more than 

50 MeV larger than 2637 MeV. Thus we conclude that, at present, we find no fully convincing understanding 
00 ■ of the quantum numbers of the observed resonance. 

ON 

Recently, DELPHI has observed a narrow resonance D(2637) in the D*tttt channel with a total width of less than 
the detector resolution: 15 MeV 0. The mass of the observed resonance turns out to be in perfect agreement 
with predictions of the quark models ^ [| for the charm radial excitation D*'(J P = 1~). This coincidence has 
■ lead to a quick identification of the discovered resonance state with the radially excited vector charm meson. 

In this letter we reconsider the identification of the quantum numbers of the observed resonance by submitting 
it to the following criteria: 

ON 



i the resonance mass should be ~ 2637±6 MeV; 

ii the resonance width should be < 15 MeV; 



iii the resonance should have a sizeable branching ratio of the channel £>(2637) — > D*-ktt in which it has been 
observed. 

We estimate the decay rate of a radial excitation of the reported mass and find that, although the partial 
widths are strongly model-dependent, the total rate is conservatively estimated to be significantly larger than 
about 50 MeV. We thus conclude that the observed width of only < 15 MeV is hardly compatible with its 
identification as a radial charm excitation. 

An identification of the observed resonance with an orbital charm excitation seems to be more favourable: 
two mesons with the quantum numbers J p = 2~, 3 _ (quark orbital momentum L = 2) could have a width of 
the observed order of magnitude. However, theoretical estimates yield a mass for the latter orbitally-excited 
states approximately > 50 MeV above the reported value. In addition, in this case, one would expect two 
neighbouring states with a mass difference of about 30-50 MeV (similar to the mass difference between D 2 + (2460) 
and Z?!+(2420) in the j=3/2 positive-parity sector), while only one resonance has been reported. Altogether 
we conclude that, at the moment, there is no fully convincing understanding of the quantum numbers of the 
observed resonance. 

Our analysis is based on combining the heavy-quark symmetry relations for the transition amplitudes between 
heavy mesons through the emission of the light hadrons with the quark-model estimates. Namely, we estimate 
the decay rates of radially and orbitally excited charm resonances into Dn, D*tt, Dtttt and D*tttt assuming the 
resonance mass of 2637 MeV as measured by DELPHI. 

To obtain estimates of the branching ratios of the three-body decays with two pions in the final state we 
treat them as cascade two-body decays .0(2637) — > (D,D*)R —> (D.D*)mr through the intermediate Breit- 
Wigner resonance with relevant quantum numbers. The status of the (7r7r)/ = o,.r=o channel is not well-defined, 
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and we varied the corresponding a resonance mass in the range 400 ~ 800 MeV and the width in the range 
r(er) ~ 700-^900 MeV. The low-energy (tttt)i = i partial-wave is dominated by p. Higher partial waves of the irir 
system give negligible contributions. 

Table [l] lists the candidate charm states and their allowed decay modes as given by the spin-parity conser- 
vation. 

For heavy-meson decays, additional constraints are given by the heavy-quark (HQ) symmetry [Q. Namely, 
in the heavy quark limit the heavy quark spin decouples from other degrees of freedom and remains conserved 
in hadron transitions. Thus, in strong decays of heavy hadrons, the total angular momentum of the heavy and 
light degrees of freedom are conserved separately, in addition to the conservation of the parity and total angular 
momentum. Hence, with respect to strong decays, heavy hadrons can be assigned an additional conserved 
quantum number, j, which is the total angular momentum of the light degrees of freedom. The consequences 
of the HQ symmetry for hadron transitions have been worked out by Isgur and Wise Namely, the HQ 
symmetry allows one to relate to each other different amplitudes of strong transitions between the states with 
fixed j and j', the latter being the angular momenta of the light degrees of freedom in the initial and final 
hadronic states, respectively. Table |l| also presents the HQ symmetry allowed transitions in terms of the few 
independent amplitudes. The 0(1/ 'ttlq) corrections in the effective Hamiltonian yield corrections to these HQ 
symmetry relations. However, for our order-of-magnitude analysis these corrections are generally unimportant 
and will ne neglected unless explicitly specified. 

For the calculation of the independent amplitudes (a, /?, 6, and £ in Table 0) one needs a non-perturbative 
approach. We apply here a naive quark-pair-creation ( 3 Po) model || which, in spite of its simplicity, has proven 
to provide a reasonable quantitative description of the two-body hadronic decays. The model is based on the 
assumption that the the spectator quarks do not change their SU(3) quantum numbers, nor their momenta 
and spins. The created quark-antiquark pair should be therefore in a 3 Pq (J pc = ++ ) SU(3) singlet state of 
zero total 3-momentum. More details concerning the model can be found e.g. in Ref. |(|. To compute the 
transition amplitudes one needs an overall transition-strength constant which determines the amplitude of the 
production of the light qq pair from the vacuum and the wave functions of the initial and final mesons. For the 
overall strength constant we use the value 7 = 2.2 as found from the analysis of the hadronic decays in the light 
sector H . For the nonperturbative meson wave functions we assume an harmonic-oscillator approximation. 
The remaining parameters to be fixed are the size of the light-light and heavy-light wave functions. They 
are given by two radii the light-light (i?) and heavy-light (Rd)- For the ground state these radii are simply 
R 2 = 2/3 < (f q - r q f >. 

The heavy meson has to be smaller than the light one. From estimates with different potentials, we assume 
the radii to satisfy the relation R 2 D / R 2 ~ 0.5-^0.7 and allowed R 2 — 6-^9 GeV 2 : such values of R are compatible 
with previous descriptions of the spectrum and decay rates Q and in addition, we check, Table ||, that the 3 Pq 
model with these parameters describes correctly the experimentally observed D 2 + — > [D, D*)ir decay rates. For 
the D* — > Dir transition, Table 0, we also express the decay width in terms of a dimensionless coupling constant 
defined as follows 

(Lfi(p2)r + {q)\D*+(p 1 ))=g D . D ^q li ^, (1) 

where e IM (pi) is the vector-meson polarization vector, the states being normalized covariantly and where we have 
omitted the momentum conservation delta function. The gD*Dn of the 3 Po model agrees with the experimental 
bound although seems to be a bit small compared with other theoretical estimates. Notice however that the 
3 Pq model is non-relativistic and as such does not describe properly the soft pion limit. In the D* — * Dir decay 
the pion is produced almost at rest and the model is indeed expected to underestimate the coupling constant 
as observed. But the model has proven to work rather well for hadronic resonance decays in the usual domain 
of the emitted pion energies, i.e. E„ ~ 300-600 MeV. Indeed, as we will see later, the model estimates for the 
decay of a heavier radially-excited D' are in better agreement with covariant methods. 

Having thus fixed the ranges of the basic parameters from the light sector and D 2 + — » (D, D*)ir decays, we 
apply the model to the analysis of the decay of radially and orbitally excited negative-parity states. 

Radial excitation D\- 

The decay of the vector radially-excited D x - into Dir and D*w is governed by the following amplitudes 

(D°(p 2 )ir+(q)\D*'+( Pl )} = gn^D^A, (2) 
(D* (p 2 )ir+(q)\D*'+( Pl )} = igo^D^^fiV^A 
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with the coupling constants gD*'Dn and gD*'D*n to be determined on the basis of a dynamical approach. In 
the heavy-quark limit the constants are related to each other as follows 

9d"Dk = Mb*' Qd* ' d* 7t ■ (3) 

To estimate these coupling constants one can applied various theoretical approaches. For instance, one may 
use the PCAC definition of the pion field (although the pion is not soft at all in this decay) in which case a 
complicated problem of calculating the coupling constants of interest is reduced to a relatively simpler one of 
calculating the D*' — > (D,D*) transition form factors through the axial-vector current. Namely, one finds 

gvpn = -j- [(My + M P )A 1 (0) + (My - M P )A 2 (0)] . (4) 

JTT 

For estimating the meson transition form factors / and a+ we used the relativistic dispersion approach of Ref. 
H adopting wave functions of the ground-state and radially-excited D which provide the values fr> — 200 MeV, 
f D ~ 240 MeV and f D *. < 400 MeV. 

Actually, one observes a suppression of the form factor A\ in the D*' — > D transition, as compared with 
the D* — > D one, because of the orthogonality of the wave functions of the orbitally-excited and the ground 
states : namely, Af '^(O) ~ 0.5 and *'^ D (0) ~ 0.2 - 0.3. In the absence of the second term of (§) the 
D* 1 — » (D,D*) transition would thus be suppressed by a factor 2 to 4 in rate. However such a suppression 
due to orthogonality in the soft pion limit is expected to be reduced when the momentum recoil (which is 
large in this case) is taken into account. The transition D*' — ► (D,D*)tt would only produce a soft pion if 
the mass of the D*' would be close to that of the D^*\ thus canceling the second term in (^). But we are 
not in such a situation, the second term is not negligible since A® ""^(O) ~ 0.9 — 2.5 and one ends up with 
the relation gw' Dtt — (0.5 — 1.5)gD* Dtt where gvD-x — 15. Notice that the large uncertainties in gD*'D-K are 
connected with a strong sensitivity of the latter to the subtle details of the sign-changing wave function of the 
radially-excited state. Altogether, and using the HQ symmetry relation (|^) we find a rather large range for our 
estimate: T(D*' -> Dtt) = 20 - 200 MeV and T(D*' -> D*ir) = 25 - 250 MeV. 

To obtain more precise estimates we also used the 3 Pq model and found values in the region gu-'D-K = 13-j-15 
which are compatible with PCAC based estimates. Table ^ presents the 3 Pq model estimates of the D' x _ (2637) 
decay rates. 

So for the sum of the decay rates of the channels T(D*' — > Dtt) + T(D*' — > D*tt) one finds rather uncertain 
estimates ranging from 40—50 MeV, which is only slightly above the reported width of the DELPHI resonance, 
to several hundreds MeV, which is far above. 

Notice that anyway an important contribution to the total width of the radially-excited D*' is given by the 
decay channel D*' — > D*(irir)i=o allowed by the HQ symmetry in the S-wave (Table Q). The 3 Pq model estimate 
of its rate is 120-160 MeV. So, conservatively, one cannot expect the total width of the radially excited state 
D *',_ (2637) to be less than 50 MeV and presumably it should be much broader. This has the double effect of 
predicting a large branching ratio for the observed D*ixix channel, but unhappily also a large total width. Thus 
we conclude that identification of the resonance D(2637) as a radial excited J p = 1~ charmed meson is hardly 
compatible with the reported value of the resonance total width < 15 MeV. 

Orbital excitations D 2 - 3- 

To proceed with these states, we assume the HQ symmetry at the level of the transition amplitudes but use 
the physical masses to compute the relevant phase-space factors. Namely, we calculate the transition amplitudes 
of the modes D 3 -j— 5 / 2 — * D^n and D 3 - j— 5 / 2 — > D(*\mr)i=o,i and determine all other related amplitudes 
through the HQ symmetry relations listed in Table [l]. 

Notice that the transition of D2-J—5/2 into the four positive parity states D2-J—5/2 ^1/2 3/2+"' — * 
7T7T, is allowed, in the heavy-quark limit, only in the D-wave, since the latter positive parity states with 
j = 1/2 or j = 3/2 cannot be produced with an S-wave pion from a j = 5/2. A rough estimate convinced us 
that, due to the small final momenta, these decays will be even more suppressed than the other D^tttt channels 
to be discussed later. 

Taking into account the phase-space factors yields the results listed in Tables || and ^. Some comments 
on the presented numbers are in order. It can be seen from Table [j] that the decay into D(*) (ttk)i=i can go 

3 We are indebted to Damir Becirevic and Alain Le yaouanc for attracting our attention on this suppression which they find 
even stronger in a related approach |10[ based on the Dirac equation. 



3 



through a P-wave between tttt and the charmed meson. Compared to the F-wave for DMtt and the D -wave for 
D^*' (tttt)i = o, this decay channel has to be considered even though only a small part of the p-meson Breit-Wigner 
tail is included in the phase space. The Z)W (7T7t); = o is suppressed by one order of magnitude as compared to 
the £>W( 7T7r); = i due to the centrifugal barrier suppression. 

The physical D 3 - state is dominantly the D 3 -j =5 / 2 one, and one can safely estimate r(D 3 - ) ~ T(D 3 - ^=5/2), 
since a possible small admixture of the narrow D 3 -j =7 / 2 state does not practically change its width. For the 
D 2 - state however this is not the case: one might expect a sizable increase of the width of the physical D 2 - 
with respect to the width of the D 2 - j =5 / 2 . For example, even a small admixture of the D 2 -j =3 / 2 to the 
dominant D 2 - j =5 / 2 may increase the total width of the physical D 2 - meson, since one expects T(D 2 - j = 3/ 2 ) 3> 
r(D 2 - ^=5/2) due to the P-wave decay mode D 2 - ,^=3/2 - * D*ir allowed by the HQ symmetry. 

A similar situation has been observed in the positive-parity Dj= 3 / 2 multiplet: namely, the experimental 
ratio of the decay rate of Di+(2420) which is dominantly -Di+ J= 3/2 and the decay rate of D2+(2460) which is 
practically the rate of a pure D 2 + j =3 / 2 is: 

r( J D 1+ (2420))/r(£)2+(2460)) = 0.71 (5) 

that is about twice larger than the leading-order HQ symmetry estimate 

r( J D 1+J=3/2 )/r(£> 2+J=3/2 ) = 0.3. (6) 

The latter value is obtained by assuming the leading-order HQ symmetry relations between the amplitudes and 
taking the physical masses of the corresponding states for the calculation of the relevant phase-space factors. 

This discrepancy can be solved by invoking the l/m c corrections e.g. by assuming a small admixture of a 
broad -Di+ J= i/2 to a narrow -Di+ J= 3/2 in the physical D!+(2420) as proposed in Ref. JTlJ . Another possibility 
(see and refs therein) is to have a rather strong increase in the decay rate of a pure Di+ j—3/2 due to the 
l/m c corrections to the effective Hamiltonian (A combination of these two variants is of course also possible). 
Anyway the Dj= 3 / 2 positive-parity sector prompts that the net effect of the l/m c corrections is a doubling of 
the ratio T(D J=j _ 1/2tJ )/T(Dj =j+1/2 j). 

Hence the obtained value of the D 2 - width should be considered as a lower bound of the physical D 2 - 
width. Still, from the comparison with the j = 3/2 sector, we expect the physical D 2 - width not to exceed the 
HQ estimate of the D 2 - j =5 / 2 state by much more than a factor 2. 

Similarly, the subleading l/m c effects can influence also the rate of the transition D 2 - ,j=s/2 ~* ^j=i/2 3/2+ 7r - 
A rough estimate based on the eqs. (||), @, where an S-wave 0(l/m c ) decay leads to an increase of about 10 
MeV of the width, and considering that the phase space is smaller in the decay of a D(2634) into positive parity 
resonances than in the decay of the latter into the ground state, we are not too worried. Still, a closer scrutiny 
of these 0(l/m c ) effects would be welcome. 

Finally, we conclude that the orbitally excited D 2 - and D 3 - charm resonances with the mass in the region 
of 2640 MeV can have the width of the order reported by DELPHI. However, an identification of the resonance 
13(2637) with an orbital excitation in the charm system is not straightforward since the reported mass seems 
to be significantly smaller than the theoretical expectations. 

For example, the Godfrey-Isgur (GI) model which describes with a good accuracy nearly all known mesons 
predicts the mass of the D 3 - to be 2830 MeV. Taking into account that for the D 2 + state the GI model gives 
2500 MeV which is 40 MeV heavier than the observed value of 2460 MeV, we can expect for the D 3 - state the 
mass < 2800 MeV. From the typical mass-splitting between the states with neighbouring j = 2 and j = 3 in 
the GI model, one could expect the D 2 - mass near 2750 MeV. 

The quark-gluon string model |l3| has predicted the masses of the charm resonances M(D 2 - ) = 2660 ± 70 
MeV and M(D 3 -) = 2760 ± 70 MeV. The partial rates of the latter were estimated to be r(£> 3 - -> Dn) = 
1.3 2.0 MeV and r(L> 3 - -> D*n) = 3.5 -i- 7.0 MeV yielding the total width of D 3 - in the necessary range. 
On the other hand a large mass splitting between the 2~ and 3~ charm states signals that the D 2 - in the 
quark-gluon string model contains a big admixture of the lighter state D 2 - j =3 / 2 - So, for the mass of the pure 
D 2 - j=5/a state one would expect a higher value. 

Altogether, from the above theoretical analyses we could expect the mass of the D 2 - j =5 / 2 in the region of 
2650 ~ 2750 MeV which is only marginally compatible with the reported resonance mass of 2637 MeV. 

The branching ratios of the D*tttt channels in table 5 and 6 are in the range of a few percent, which may 
seem a little small for these channels to have been observed. Still some non-resonant irn contributions, which 
are difficult to estimate, will add up. 



4 



Summing up, if the resonance £>(2637) with the width of < 15 MeV is confirmed by further analyses (at the 
moment CLEO and OPAL do not see it |14[]) then the theoretical understanding of its quantum numbers is not 
clear: in spite of the coincidence of the observed mass with the predicted mass of the radial charm excitation 
D' x _, its interpretation as a radial excitation is completely ruled out by the small observed width. 

On the other hand, although an identification of the state with the D 2 - orbital excitation seems to be 
appropriate from the viewpoint of the total width, its mass seems to be too low compared with the quark-model 
theoretical estimates. Its D*mr branching ratio of a few percent is a little small. More importantly in the case 
of the orbital excitations in this mass region, one would expect two neighbouring resonances with the width of 
order several MeV each and with the mass difference of about 30—50 MeV, corresponding to the D 2 - and In- 
states. The published plots do not show any sign of a neighbouring 3~ resonance. 

Thus we conclude that the experimental confirmation of this resonance would put forward a challenge of its 
proper theoretical understanding, unless a neighbouring slightly heavier resonance was found. 
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Table 1: Decay modes of possible candidate states allowed by spin-parity conservation. Modes allowed also 
by the HQ symmetry, i.e. corresponding also to a separate conservation of the total angular momentum of the 
light degrees of freedom, are listed in bold and HQ symmetry relations for the corresponding amplitude squared 
from la] are given (without phase-space factors included). 







D 2 -,j=5/2 


D 3-,j=5/2 


Dtt 
D*ir 


L=l[f/3 2 ] 


L=l,3[a 2 ] 


L=3[fa 2 ] 
L=3[fa 2 ] 


D**ir: 

D 0+,j = l/2^ 
D\+ J = l/2 7r 
D l + J=3/2 7r 
D 2+,j=3/2^ 


L=0,2 
L=0,2 
L=2 


L=2 
L=2 
L=2 
L=0,2,4 


L=2,4 
L=2,4 
L=2,4 


£>*(7T7r) ;=0 


L=0,2 


L=2[§<5 2 ] 
L=2[|<5 2 ] 


L=2[<5 2 ],4 


D{-ktt)i =1 
D*(itit) 1=1 


L=l 

L=l,3 


L=l[f£ 2 ],3 
L=l[Ie 2 ],3 


L=3 
L=1K 2 ],3,5 



Table 2: The D* -> Dtt transition. 





Exp. 


3 Po model 


Other estimates 
(see refs in ||) 


9D'Dtt 

T{D* -> Dtt) 


< 21 

< 89 KeV 


7±1 

7H-10 KeV 


7+21 

7-^90 KeV 



Table 3: Decay rate of the D 2 + — > (D, D*)n transition. 





Exp. |' 


n 


a Pq model 


r(D7r)/r(L>*7r) 


2.3±0.9 

23±5 MeV 


2.6 

11-22 MeV 



Table 4: Decay rates of the D[_ — > (D, D*)tt transition in the 3 Pq model. 



D'(l-) 


-> Dtt 


150- 


-220 MeV 


D'(l-) 


-> L>*7T 


200- 


-300 MeV 


D'(l-) 


-> .0*7777 


120- 


-160 MeV 



Table 5: Branching ratios of the D%- ,3=5/2 decays in the 3 Pq model. 



Br(Dir) 
Br{D*iz) 


0.65-0.7 


Br{D(irir) l=0 ) 

Br{D(TTTT) l = 1 ) 


< 0.01 
0.2-0.3 


Br(D*(TTTr) l=0 ) 

Br{D*(TTTT) l = l ) 


< 0.003 
0.02-0.03 


Ttot 


6-14 MeV 



Table 6: Branching ratios of the D 3 - ,^=5/2 in the 3 Pq model. 



Br(Dir) 
Br(D*ir) 


0.65-0.74 
0.23-0.27 


Br(D(7T7T) l=0 ) 

Br{D{^)i =l ) 


<0.001 


Br(D*(TTTT) l=a ) 
Br{D*(TTTT) l = l ) 


< 0.003 
0.03-0.06 




8-22 MeV 



6 



